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We report the measurement of the anisotropic AC polarizability of ultracold polar 40 K 87 Rb 
molecules in the ground and first rotationally excited states. Theoretical analysis of the polar- 
izability agrees well with experimental findings. Although the polarizability can vary by more than 
30%, a "magic" angle between the laser polarization and the quantization axis is found where the 
polarizability of the \N = 0,m^ = 0) and the \N = l,mjv = 0) states match. At this angle, rota- 
tional decoherence due to the mismatch in trapping potentials is eliminated, and we observe a sharp 
increase in the coherence time. This paves the way for precise spectroscopic measurements and 
coherent manipulations of rotational states as a tool in the creation and probing of novel quantum 
many-body states of polar molecules. 



PACS numbers: 03.75.-b, 37.10.Pq, 67.85.-d, 33.20.Bx 

The creation of a gas of ultracold polar molecules with 
a high phase space density pQ brings new possibilities 
beyond experiments with ultracold atomic gases. In 
particular, long-range, anisotropic, and tunable dipole- 
dipole interactions open the way for novel quantum gases, 
with applications including strongly correlated many- 
body systems [2HS], precision measurement [6 , and ul- 
tracold chemistry [7J [8]. Molecules also have complex 
internal structure with many more internal degrees of 
freedom than atoms. In particular, the rotational degree 
of freedom provides a set of long-lived excited states that 
are easily coupled to the ground state with microwaves. 
Because of the accessible frequency in the microwave do- 
main and the narrow intrinsic linewidth, the transition 
between rotational states could be ideal for use as a spec- 
troscopic probe of the system. For example, such a nar- 
row transition could be used to measure small energy 
shifts due to dipolar interactions [9]. In addition, dipole- 
dipole interactions can be realized without applying a 
DC electric field but instead by directly coupling the two 
lowest rotational states with a microwave field. Within 
the rotating frame of the microwave transition, there is 
a strong dipolar interaction, which can be used to model 
novel quantum many-body Hamiltonians [10] • Addition- 
ally, the microwave detuning and power can be varied to 
modify collision dynamics [11 . This has been proposed 
as a way to achieve a topological superfluid of paired 
fermionic polar molecules [T2] . 

A prerequisite for such experiments is long coherence 
times for the interaction between a microwave field and 
the rotational states. However, for molecules confined in 
an optical dipole trap, the difference in AC (or dynamic) 
polarizability between different rotational states must be 
considered [13] . A difference in polarizability leads to 
different trap frequencies and spatially dependent vari- 
ations in the rotational transition frequency, which can 
lead to dephasing and decoherence [14]. For atoms, if 
the trapping light is far detuned compared to the en- 



ergy splitting between two states, their AC polarizabili- 
ties will be nearly equal. In contrast, for molecules, the 
ground and the first rotationally excited state with the 
same angular momentum projection onto the quantiza- 
tion axis have different parity and will therefore couple 
to different electronic excited states. This can result in 
AC polarizabilities that differ by more than 30%, even 
when the light is far detuned [15] . 

In atomic systems, it is possible to adjust the wave- 
length of the trapping light such that the polarizabili- 
ties of two states of interest (often clock states in alka- 
line earth atoms) are the same [13] . Although, in prin- 
ciple, one could find such a "magic" wavelength trap 
for molecules [16 , the large number of additional states 
from rotation and vibration makes it difficult to find 
a suitable wavelength that is sufficiently detuned such 
that off-resonant light scattering is negligible. However, 
molecules provide a different way to adjust the polariz- 
ability; the AC polarizability of a molecule depends on 
the relative orientation of the molecule and the polariza- 
tion of the trapping light [T5] . 

In this Letter, we explore the interaction between the 
trapping light and the molecules by examining the real 
part of the AC polarizability. (The imaginary part of 
the polarizability was reported in a previous study of 
the lifetime of molecules trapped in a three-dimensional 
lattice [17].) In particular, we determine how the po- 
larizability of the rotationally excited states depends 
on the relative orientation of the quantization axis, z 
(which in these experiments is given by a magnetic field 
of 545.9 G), and the polarization of the trapping light. 
A "magic" angle exists where the polarizability of the 
\N = 0, 77i 7v = 0) (where N is the rotation quantum 
number, and tun is its projection onto the quantiza- 
tion axis) and \N = l,rn/v = 0) states match, mak- 
ing the AC Stark shift the same for these two internal 
states of the molecule. This "magic" angle is expected 
at cos 2 = 1/3, or ~ 54 degrees [18 . We observe 
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FIG. 1: (a) Experimental schematic. The lattice beam prop- 
agates along x , the magnetic field points in the z direction, 
and the polarization of the lattice light makes an angle with 
the magnetic field in the y-z plane, (b) Schematic of rota- 
tional energy states. The degeneracy of the N — 1 level is 
split in a magnetic field, (c) A sketch of the optical dipole 
potentials for the |0, 0) and |1,0) states. A Gaussian is over- 
laid to show the density distribution of the molecular cloud 
in the trap. When the two states are connected by a 2.22 
GHz microwave drive, there is effectively a spatially varying 
detuning across the cloud due to the difference in the trap 
potentials. 



state is broken by the hyperfine interaction, specifically, 
the interaction between the nuclear quadrupole moment 
and the rotation of the molecule [I9j [23] . In the rotation- 
ally excited states, this coupling gives a specific nuclear 
spin state up to a 5% admixture of other hyperfine states; 
however, we will ignore these contributions and will work 
only with states whose dominant nuclear spin character 
is mf = — 4, mf h = 1/2 unless otherwise noted, where 
mf and m^ b are the nuclear spin projections onto the z 
axis. For more discussion of the nuclear quadrupole cou- 
plings, see the online supplementary material. At a field 
of 545.9 G and with no light, we find that the |1, 1) state 
is 58 kHz above the |1, — 1) state and the |1, 0) state is 268 
kHz above the |1, 1) state (FigjlJ)). The splitting between 
the tun = ± 1 states is of the same order of magnitude 
as the AC Stark shifts induced by the trapping lasers. 
Consequently, we calculate the anglar dependence of the 
dynamic polarizability of the three rotationally excited 
states using perturbation theory that includes the light- 
induced couplings between the bare states (no light). 

Using the general formalism of Ref. [24], the complex 
dynamic polarizability for the |0,0) state is given by 



«|0,0) = oHl + 2Q ^) 



(1) 



this magic angle both through direct measurements of 
the AC polarizability and through measurements of the 
coherence time for driving the |0,0) to |1,0) transition. 
In the theory component of this Letter, we extend the 
ideas of mixing of rotational levels due to static electric 
fields of Ref. [15 to include mixing due to the intrinsic 
nuclear electric-quadrupole and Zeeman interactions dis- 
cussed in Ref. [19 . We compare our experimental results 
to the analytic results of an approximate Hamiltonian, 
which mixes the three projections of the first rotation- 
ally excited state. 

Molecules have a definite orientation that strongly af- 
fects their polarizability. In the frame of a diatomic 
molecule such as KRb, the DC (or far-off-resonance 
AC) polarizability along the internuclear axis is much 
stronger than the polarizability perpendicular to the in- 
ternuclear axis. This has been successfully exploited to 
align molecules with intense laser pulses in order to study 
the stereodynamics of chemical reactions [20-22 . The 
N = rotational ground state is spherically symmetric 
and therefore its AC polarizability has no dependence on 
the relative orientation of the laser polarization and the 
quantization axis. However, the rotational wavefunctions 
of the three projections (rajv) of the N = 1 rotation- 
ally excited state have well-defined orientations relative 
to the quantization axis. For example, the |1,0) state 
corresponds to a p z orbital that is aligned along the z- 
axis, and therefore the polarizability is the largest when 
the polarization of the AC field is along the z-axis. 

The degeneracy of the three projections of the N = 1 



where the "reduced" polarizabilities ay and a±_ are the 
parallel and perpendicular (with respect to the inter- 
molecular axis) polarizabilities that describe the aver- 
aged contributions from ro-vibrational states of all elec- 
tronically excited 1 S + and 1 Ii potentials. As mentioned 
above, the polarizability of the |0, 0) state is independent 
of the angle between z and the polarization of the ID 
optical lattice used to trap the molecules (see Fig. |TJl) . 

For the |1,0) and |1,=L1) states, the dressed or mixed 
polarizabilities at laser intensity / are given by the to- 
tal Stark shift divided by the total intensity, ctj = 
_ E j( I )- E j(Q) ^ wnere with j = 1,2,3, are the eigen- 
values of the 3x3 Hamiltonian 
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and ei, 62, and 63 are the energies for states |1,0), |1, 1) 
and |1, — 1), respectively, at 1 = 0. 
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FIG. 2: Parametric heating resonances in the far-off- 
resonance optical dipole trap for = 57 degrees. The y-axis 
shows the rms size in x of an expanded gas of KRb (Rb) after 
5 ms (21 ms) of time of flight. The curves have been offset 
vertically for clarity. Using Gaussian fits (lines), we determine 
the center of the parametric heating resonances for (from bot- 
tom to top) Rb, and KRb in the |0,0), |1,0), |1,1), and |1,-1> 
states. The resonant frequency allows us to extract the trap 
depth for each state. 
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FIG. 3: The AC polarizability of KRb at 1064 nm for the 
|0,0) (black squares), 1 1 , 0) (blue circles), |1,1) (red inverted 
triangles), and |1, — 1) (green triangles) states. Error bars are 
from the fit uncertainty in the center of the parametric heating 
resonances and correspond to ±1 standard deviation. Theory 
lines are a simultaneous fit to Eqn. [I] and the polarizabilities 
from the solution of H with three free parameters m , a±, 
and a || . Open circles represent a separate measurement where 
the polarizability is extracted from the shift in the microwave 
transition frequency. 



We measure the AC polarizability of the molecules in a 
one-dimensional optical lattice with a peak intensity, 2o> 
of 2.3 kW/cm 2 and a wavelength of A = 1064 nm [171125]. 
Using a single microwave pulse, we can selectively trans- 
fer population from the rotational ground state, |0, 0), to 
any of the projections of the N = 1 state with near 100% 
efficiency [19 j. We measure the trap depth, /TxRb? from 
which we can extract the AC polarizability, by measuring 
the parametric heating resonance (see Fig. [2|. We mod- 
ulate the intensity of the optical lattice for 4 ms with 
an amplitude of approximately 10% of the total depth. 
The modulation frequency is varied to find the resonant 
frequency where molecules are excited from the lowest 
band of the lattice to the second excited band. In the 
deep lattice limit, this resonant frequency is twice the 
trap frequency. However, we operate our lattice in an 
intermediate intensity regime where the relationship be- 
tween the trap depth and the resonant frequency must be 
extracted from a numerical solution to the lattice poten- 
tial. We determine the polarizability, without having to 
characterize the optical beam parameters such as power 
or beam waist, by comparing the molecular results to a 
similar measurement of the trap depth for Rb. The polar- 
izability of KRb is then given by a K Rb = «Rb^KRb/^Rb, 
where a Rh /h = 3.242 x 10~ 5 MHz/(W/cm 2 ) at 1064 nm 
[26] , and h is Planck's constant. 

In Fig. [3J we show the polarizabilities of the N = 
ground state as well as the three projections of the N = 1 
rotationally excited state as a function of which is var- 
ied by adjusting a half-waveplate in the lattice beam. 
Although the rotation of the half-waveplate allows us to 



choose with precision better than one degree, the ab- 
solute alignment relative to the magnetic field has an es- 
timated systematic uncertainty of ±3 degrees. We fit 
Eqn. [T] and the polarizabilities from the eigenenergies 
of H to the experimental data with three free parame- 
ters, m , an and a±. From the best fit, we determine 
that the "magic" angle m = 48(4) degrees, a\\/h = 
10.0(3) x 10- 5 MHz/(W/cm 2 ), and a±/h = 3.3(1) x 10" 5 
MHz/(W/cm 2 ). 

Hyperfine couplings between the |1,0) and |1,±1) 
states result in a small change of the predicted magic an- 
gle from 54 degrees. For the polarizabilities and intensity 
given above, we expect m = 52 degrees, which agrees 
with our measurement to within the error. For compari- 
son to the measurement, theoretical values of the polar- 
izabilities a || and a± are obtained using non-relativistic 
potentials and dipole moments described in Ref. [15] , and 
are h x 12.2 x 10~ 5 MHz/(W/cm 2 ) and h x 2.01 x 10~ 5 
MHz/(W/cm 2 ), respectively. In principle, an and a± de- 
pend on the ro- vibrational state of the molecule, however, 
we find in the calculation that for small N, the polariz- 
abilities a\\ and a± are independent of N to better than 
0.01%. 

Because the bare states are mixed by the lattice light, 
the dresssed state polarizabilities depend on the inten- 
sity of the light. This is effect is strongest for the |1, ±1) 
states, where the energy splitting of the bare states is rel- 
atively small. However, we experimentally verify that at 
the intensity used here this effect is small. For example, 
for a 50% increase in intensity with = 93 degrees, we 
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see only a 3.9(8)% decrease in the polarizability of the 
|1, 1) state. 

We can also use an independent measurement of the 
shift in the microwave transition frequency to determine 
the polarizability as a function of 0. The difference in 
polarizability between the |0,0) and |1,0) states results 
in a shift in the microwave transition frequency: 

f = fo + (a m - a m )I /h +^-^+Af, (2) 

where / is the measured transition frequency, fo is the 
transition frequency from |0, 0) to |1, 0) with 1 = (mea- 
sured after release from the lattice), Cjj/(4tt) = y^f^Sk 
is the trap zero point energy in the lattice divided by /i, 
m is the mass, Af = a|1 ' 0) ~ a|0 ' 0) k B T/h is the shift in 
the center of the transition frequency caused by the spa- 
tially dependent detuning, ks is Boltzmann's constant, 
and T = 400 nK is the temperature. The transition fre- 
quency is only sensitive to the polarizability difference 
between the |0,0) and |1,0) states, but can be compared 
with the polarizability from the direct measurement of 
the trap depth by fixing a| ,o) an d then solving Eqn.[2]for 
a|i 5 o). The resultant polarizabilities are shown in open 
symbols in Fig. [3] We see a good agreement between the 
two methods. 

To study the effect of the magic angle, we measure 
the rotational excitation coherence time as a function of 
angle, and the results are shown in Fig. [4j We measured 
the coherence time between the \N = 0,mjv = 0, mf = 
-4,m^ b = 1/2) and \N = l,m N = 0,mf -3,mf h = 
1/2) states with Ramsey spectroscopy. The probe pulse 
is 40 /is long, and the probe frequency is detuned from the 
resonance by 3 to 12 kHz. We fit the Ramsey oscillation 
as a function of time to a damped sine wave to extract 
the coherence time (see inset of Fig. |4|. 

The coherence time due to the mismatch in polariz- 
ability should scale as one over the difference in dE/dl, 
which, because of the intensity dependence of the polariz- 
ability, is not the same as the difference in polarizability. 
Technically, the Af term in Eqn. [2] should also use the 
"local polarizability", dE/dl, instead of a because the 
trapped gas experiences only a small range of intensities. 
However, in Eqn. [2] the difference between a and dE/dl 
gives a small correction to a term that accounts for less 
than 10% of the total frequency shift, and is therefore 
negligible. On the other hand, because the coherence 
time depends critically on this polarizability difference, 
the "local polarizability" must be used. 

We fit to the data with a simple model that includes 
the angular dependence: 
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FIG. 4: The Ramsey coherence time measured in the one- 
dimensional optical lattice as a function of angle. A sharp 
increase in coherence time is observed at the "magic" angle 
where the polarizabilities of the |0,0) and 1 1 , 0) states are 
matched. Inset: A Ramsey oscillation fit to a damped sine 
wave to extract the coherence time for = 51 degrees. 



spread of energy across the cloud in the |0,0) state, and 
\dE\ ^/dI — dE\i^/dI\ depends upon ay, a±, /, ei, 
62, and 63. Using the measured values of /, ei, 62, and 63 
and the fitted values of ay and a± from Fig. [3| we use 
Eqn. [|and obtain T 2 = 1.5(2) ms, AE/h = 3.3(4) x 10 4 
s _1 , and m = 46.5(5) degrees. The expected value of m 
for our values of /, a±, and a\\ is 48 degrees, which agrees 
with the fitted value to within the systematic error. Note 
that the expected value of m is different when probing 
the difference in dE/dl rather than a. Possible sources 
of T2 range from technical noise to resonant dipole-dipole 
interactions [9], and further study is required to under- 
stand this limit to the coherence time. 

In conclusion, we have measured the angular depen- 
dence of the AC polarizability of ultracold KRb and ob- 
served the "magic" angle where the polarizability of the 
|0,0) and |1,0) states match. At this angle, we are able 
to increase the coherence time between the two rotational 
states by an order of magnitude. This opens the way for 
the use of the rotational states for precision spectroscopy 
as well as the modification of the collision dynamics in 
optically trapped molecular samples. 

We acknowledge financial support from NIST, NSF, 
DOE, AFOSR, and DARPA. S. A. M. acknowledges 
funding from the NDSEG. We thank K. R. A. Hazzard, 
A. M. Rey, J. L. Bohn, and G. Quemener for useful dis- 
cussions. 



T = IU (1/T 2 )2 + 



dE im /di-dE M /di 

dE m /dI 



AE/h) , 



(3) 

where T2 is the coherence time from all sources of 
decoherence other than the polarizability, AE is the 
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Supplementary Information 

As was shown in Ref. pQ, the nuclear quadrupole mo- 
ment couples to the rotation of the molecule and mixes 
rotationally excited states with different hyperfine char- 
acter. It is precisely this coupling that breaks the degen- 
eracy of the three projections of the N = 1 state. These 
energy shifts move the magic angle from the zeroth order 
value of 54 degrees to 52 degrees. However in the main 
text, other than to include this broken degeneracy we 
completely ignored the couplings between different hy- 
perfine states in the perturbation theory calcuations. To 
justify this assumption we now compare the results from 
the approximate 3x3 Hamiltonian in the main text to 
an "exact" model that includes the nuclear quadrupole 
couplings. 

We extend the ideas of mixing of rotational levels due 
to a static external electric field of Ref. [2 to include 
mixing of rotational-hyperfine levels due to the intrin- 
sic nuclear electric-quadrupole interactions introduced in 
Ref. pQ. First, we construct the full Hamiltonian for the 
rotational-hyperfine levels labeled by \N : mjsr, mf, mf h ), 
where N is the rotation quantum number, tun is its 
projection onto the quantization axis, mf is the pro- 
jection of the K nuclear spin onto the magnetic-field 
axis, and m^ b is the projection of the Rb nuclear spin 
onto the magnetic-field axis. It includes the nuclear 
Zeeman interaction, — g a /iAr/ a £? for atom a = K or 
Rb, where g a is nuclear g- factor and hn is the nu- 
clear magneton of atom a. It also includes the nuclear 
quadrupole interaction, proportional to {eqQ) a / {I a {I a — 
1)) E m (- 1 ) m ^2m(^, 0)T 2 ,_ m (/ a , J a ), with coupling con- 
stants (eqQ) a for each atom a that couples its nuclear 
spin to rotational states. Here, Ci m (6,(j)) is a spherical 
harmonic and T2 m (I a , I a ) is a rank-2 tensor created from 
the spin J , e is the proton charge, q is the electric field 
gradient, and Q the nuclear quadrupole moment. Fi- 
nally, we include a polarizability interaction Hamiltonian, 
— (a||0|| -\-a±0±)I, with strengths a\\ and a±, tensor op- 
erators On and Oj_ that depend on light polarization and 

rotational angular momentum N, and the laser intensity 
of the trapping light /. 

This Hamiltonian couples (1 + 3)x9x4 = 144 chan- 
nels | TV, ttitv, TRj, m^ 3 ) and has four parameters: the 
quadrupole interaction constants for each of the two 
atoms and the "reduced" polarizabilities a\\ and a±, 
which are the vibrationally- averaged v = parallel and 
perpendicular polarizabilities that include contributions 
from all excited 1 E + and x n potentials, respectively We 
find transition energies by diagonalizing the Hamiltonian 
and analyzing its eigenfunctions to connect to the states 
that have been observed experimentally Eigenstates are 
identified by the channel state with the largest contribu- 
tion. Typically this contribution is more than 90% and 
we are justified in labeling eigenstates with only one par- 
ticular set of nuclear spin projections \N,rriN, mf, mf h ). 
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TABLE I: The experimentally measured rotational-hyperfine 
transition frequencies of the lowest vibrational level of the 
X*£ + potential of KRb at zero laser intensity (i.e. without 
trapping light) and a bias magnetic field with strength B = 
545.9 G. Transitions start at the \N — 0,rajv = 0,mf = 
1/2) state and go to three hyperfine states \j) 
within the N = 1 manifold. 



Ii> 


Frequency (kHz) 


|1,0, -4,1/2} 
11,1,-4,1/2} 
11,-1,-4,1/2} 


2 228 110(1) 
2 227 842(1) 
2 227 784(1) 
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FIG. 5: The angle-dependent polarizability for the "exact" 
model (dashed line) and the approximate 3 Hamiltonian from 
the main text. 



ing new measurements of the transition energies for three 
hyperfine levels of the N = 1 rotational state where the 
effects of the trapping potential have been removed by 
measuring the transition energy directly after turning off 
the trapping light. These transition energies are given in 
Table [| 

We optimize the quadrupole interaction constants to 
fit to the experimental energies of Table [l| We find that 
£ e //i=l. 1139514(5) GHz, (eqQ) K /h=0A52(9) MHz, and 
(eqQ) Rh /h = -1.308(9) MHz. Here, B e is the rotational 
constant and h is Planck's constant. The first two val- 
ues are in good agreement with the previous estimate of 
Ref. pQ. The (eqQ)^ coupling constant has changed by 

The polarizability of eigenstate j with energy Ej(I) 
is defined as the total Stark shift divided by the total 
intensity, aj = — Ej ^~ E ^^ . For comparison to the ap- 
proximate 3x3 Hamiltonian used in the main text, we 
show the results of the "exact" model using the fitted val- 
ues of ay and a±. Figure [5] shows the polarizability as a 
function of based on theP'exact" Hamiltonian (dashed 
lines), and the approximate Hamiltonian (solid lines) for 
the four states studied here and laser intensity 7=2.35 
kW/cm 2 . From this comparison, we conclude that the 
difference between the two models is sufficiently small 
that the use of the simplified approximate Hamiltonian 
is justified given the current experimental precision. 



The two quadrupole interaction parameters were esti- 
mated in Ref. [1 based on measurements of transition 
frequencies between sub-levels of the N = and N = 1 
states in the optical dipole trap. However, when these 
measurements were made, the state-dependent polariz- 
ability and the effects of the trapping potential on the 
transition frequency were not yet understood, and the 
measurements for the different states were not taken at a 
constant intensity. Here we improve these constants us- 
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